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ABSTRACT 



Aims. Observations of high-energy y-rays recently revealed a persistent source in spatial coincidence with the binary system 77 Carinae. 

Since modulation of the observed y-ray flux on orbital time scales has not been reported so far, an unambiguous identification 

was hitherto not possible. Particularly the observations made by the Fermi Large Area Telescope (LAT) posed additional questions 

regarding the actual emission scenario. Analyses show two energetically distinct components in the y-ray spectrum, which are best 

described by an exponentially cutoff power-law function (CPL) at energies below 10 GeV and a power-law (PL) component dominant 

at higher energies. 

Methods. The increased exposure in conjunction with the improved instrumental response functions of the LAT now allow us to 

perform a more detailed investigation of location, spectral shape, and flux time history of the observed y-ray emission. 

Results. We detect a weak but regular flux decrease over time. This can be understood and interpreted in a colliding-wind binary 

scenario for orbital modulation of the y-ray emission. We find that the spectral shape of the y-ray signal agrees with a single emitting 

particle population in combination with significant absorption by y-y pair production. 

Conclusions. We are able to report on the first unambiguous detection of GeV y-ray emission from a colliding-wind massive star 

binary. Studying the correlation of the flux decrease with the orbital separation of the binary components allows us to predict the 

behaviour up to the next periastron passage in 2014. 

Key words. Gamma rays: stars - Binaries: general - Stars: binaries 



1. Introduction 

Owing to its dramatic history (the 'great eruption' of 1843, sub- 
sequent disappearance, envelopment in a huge bipolar nebula, 
etc.), 7] Carinae (77 Car) has been of interest for astronomical and 
astrophysical rese arch for more tha n a century. Rece nt observa- 
tions b y AGILE dTavani et alJ2009h and Fermi-LAT (lAbdo et al.1 
l2010bl) established high-energy y-ray emission in spatial coin- 
cidence with the location of 77 Car. Physical association, e.g. 
through observation of correlated multiband variability or orbital 
modulation of the y-ray signal, has not yet been confidently pre- 
sented. If this identification can be established, the 77 Car sys- 
tem would be the first massive star colliding-wind binary system 
(CWB) in the y-ray sky, in contrast to known y-ray binaries, 
which all contain a compact object. 

It is widely accepted (e.g., Damineli et al .1120081) that 77 Car 
is a binary system consisting of two massive stars (one being 
a member of the rare stellar class of luminous blue variables 
(LBV), the other being an O or WR star) with an orbital period 
of 5.54 yr. Stellar and orbital parameters as determined by vari- 
ous studies are listed in Table [T] Both massive stars in the 7/ Car 
system are expected to produce powerful stellar winds. Mass- 
loss rates and terminal wind velocities (as given in the table) 
are thought to be sufficiently high to form a wind-wind colli- 
sion zone of shoc ked, hot gas, wherein p article accelerati on (as 
illustrated by e.g.. lEichler & Usovlll993HDoughertv et al.ll2003h 
in general and s ubsequent y-ray emission in particular (as de- 
scribed by e.g., iReimer et al"] 120061; iBednarek & Pabichl I20TH) 
can occur. The conditions in the wind-wind collision zone de- 



pend on the orbital phase of the binary system. Therefore, these 
models predict y-ray emission that is generally modulated on 
orbital time scales. In addition there may exist complex depen- 
dencies, such as spectral cutoffs in the underlying particle spec- 
trum due to excessive energy losses, anisotropy in the Inverse 
Compton process, or types of particle transport. 



Previous analyses of Fermi-LXY data (Ab do et all 12010b; 
Far mer et"aD 1201 ll) have ela borated on the cat alogued Fermi- 
source 1FGL J1045.2-5942 (I Abdo et alJl2010ah that was found 
to be spatially consistent with the location of 77 Car. Following 
the nomenclature of the new Fermi-LAI two-year catalogue 
dNolan et al.ll2012h . the source designation 2FGL J1045.0-5941 
is used from here on. The high-energy spectrum has been best 
described by an exponential cutoff power-law function (CPL) 
that appears to be dominant below 10 GeV, and an additional 
component described by a power-law function (PL) at higher 
energies. Neither study presented a flux time history that would 
unambiguously relate the observed y-ray emission to a phys- 
ical origin in the 77 Car system. Whereas there is evidence 
for nonther mal emission of 77 Car at energies well be low the 
MeV range (iLevder et al.l 120081: ISekiguchi etail l2009h , no y- 
ray emission at energies above the energy range of the Fermi- 
LAT instrument has been report ed from the vicinity of 77 Car 
(IHESS Collaboration et alJl2012b . 
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Parameter 


Value 


Reference 


d 
P 

i 
e 

O 

a 

M A 
M B 

M A 
M B 


2.3+0.1 kpc 
2024+2 d 

45° 
0.9 

27° 
0-30° 

15.4 AU 
90 M 
30 M 

10- 3 Moyr 1 

2.5 x 10" 4 Moyr 1 

10- 5 Moyr 1 

1.5 x 10- 5 Moyr 1 

500 kms" 1 
3000 kms" 1 


Davidson & Humphreys (1997) 
Corcoran et al. (2005) 
Okazaki et al. (2008) 
Smith et al. (2004) 

Okazaki et al. (2008) 


Parkin et al. (2009) 

Corcoran (2001) 
Hillier etal. (2001) 
Verner et al. (2005) 

Hillier etal. (2001) 
Pittard & Corcoran (2002) 

Pittard & Corcoran (2002) 


Parkin et al. (2009) 


Voo,B 


Hillier etal. (2001) 
Pittard & Corcoran (2002) 



Table 1. Several important stellar, orbital and stellar-wind pa- 
rameters of the 7] Car binary system and their values as deter- 
mined by various studies. P and d are the orbital period and the 
distance from the Earth. /, e and a are inclination angle, eccen- 
tricity and semi-major axis. The parameter O represents the an- 
gle between the projection of the orbital plane along the line of 
sight and the apastron side of the semi-major axis in the prograde 
direction: <D = 270° - cd, with a> being the commonly used peri- 
astron longitude angle. M, M and Voo are stellar mass, mass-loss 
rate and terminal velocity of the stellar wind for the stars A and 
B in the system. 



2. Observations with the Large Area Telescope 

2.1. Dataset 

The Fermi-LAI started regular survey-mode observations on 
2008 August 4, providing a complete, regular sky-coverage with 
approximately uniform exposure. Accordingly, the source of in- 
terest was monitored regularly over the whole time range of the 
dataset that is used in this analysis (~35 months, 2008 August 
4 to 2011 July 5). The data were reduced and analysed using 
the Fermi Science Tools v9r23 packageQ. To reject atmospheric 
y-rays from the Earth's limb, we excluded time periods in which 
the region was observed at a zenith angle greater than 100° and 
for observatory rocking angles greater than 52°. All remaining 
photons with energy E > 200 MeV within a square of 21° base 
length (aligned in directions of right ascension a and declination 
6) centred on the nominal location of r\ Car, (a, c))=(161.265°, 
-59.685°) were used. 



2.2. Likelihood analysis 

Because r\ Car is located at low Galactic latitude and within the 
projection of the Carina- Sagittarius arm of the Milky Way, care- 
ful modelling of the diffuse y-ray emission is required. Due to 
the size of the Fermi-LAT's point-spread function (PSF) - the 
68% containment radius is more than 2° at 200 MeV - pos- 
sible contributions from nearby sources have to be taken into 
account as well. We performed maximum likelihood analyses 
using the instrument response function P7SOURCE_V6 in con- 
junction with the Galactic diffuse model galJ2yearp7v6.vO.fits 
and an isotropic background component iso-p7v6source.txt. The 
applied source model includes 38 point sources that are located 



within a radius of 15° around the source of inte rest and are 
listed in the Fermi-LAI two-year (2FGL) catalogue (JNolan et al.l 
201 ^). For each source the same spectral model as in the cata- 
logue was used. The spectral parameters of all sources within 
a radius of 10° were allowed to vary in the likelihood analysis. 
Most sources are modelled either by a simple PL 



dN A 7 ~ 



-r 



(1) 



(where A is the normalization parameter and T is the spectral 
index) or by a LogParabola function, which provides additional 
information on the spectral curvature of some sources by intro- 
ducing a break energy Eb and the parameters a and ft, 



dN _ 4(JL )-( a+ £ lo s(^>> 

dE Eb 



(2) 



Exceptions are y-ray pulsars, which are modelled by a CPL func- 
tional expression, 



dN 4 r, T , E 

— = A £" r exp( 

dE ^cutoff 



), 



(3) 



where Zs C utoff is the cutoff energy (other parameters defined as 



1 See the Fermi Science Support Centre (FSSC) website for details: 
http : //£ ermi . gsf c . nasa . gov/ssc/ 



above.) In all three cases ^ gives the differential flux. 

Source detection significance can be described by the likeli- 
hood test statistic value TS = -2ln(L maXr o/ L max ,\), which com- 
pares the ratio of two values that are obtained by a maximum 
likelihood procedure. L max $ is the likelihood for a model without 
an additional source at a specified location (the null-hypothesis) 
and L maXj \ is the maximum likelihood value for a model includ- 
ing an additional source. Surrounding sources and background 
are taken into account in both cases. The same notion can also 
be applied in comparing the likelihood of two source models 
that are described by different spectral models for a specific 
source of interest. Because the TS -value is expected to follow 
a ^-distribution for the differenc e in the numb er of degrees of 
freedom between the two models (Ma ttox et al.ll 19961) , it can be 
converted into the detection significance cr for a point source. 



3. Analysis details and results 

3.1. Spatial analysis 

In the two Galactic coordinate grids in Figure Q] the TS-values 
indicate detection significances for a trial additional point source 
(a PL with two free parameters) considered sequentially at each 
of the grid positions. The left figure covers the region around 
2FGL J1045.0-5941 at energies 0.2 to 10 GeV, the right figure 
at energies 10 to 300 GeV - each representing 35 months of data. 
The energy band has been divided according to the previous in- 
dication of two different emission components. Both TS-maps 
were obtained with the tool gttsmap. In addition, confidence con- 
tours mark the 68.3%, 95.4% and 99.7% uncertainty regions for 
the location of the source defined by the maximum TS value in 
the grids. For both energy bands the nominal position of 77 Car is 
found within the 95.4% error region. For the low band it is con- 
tained within the 68.3% error region. The formal significance for 
the low- and high-energy band is 46 cr and 12 <x, respectively. 

In the lower-energy band the y-ray signal is centred at (a, 
£)=(161.243°, -59.681°) with a fairly circular 95% confidence 
region of radius r=0.023°. In the higher-energy band the sig- 
nal is centred at (a, £)=(161.209°, -59.710°) with an elliptical 
95.4% confidence region that has a semimajor axis of a=0.053° 
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(NE-SW oriented) and a semiminor axis of Z?=0.045°. The in- 
terplay between a vast photon number at lower energies and 
statistical limitations at higher energies yields a smaller confi- 
dence region for the hard band. The angular separations d of the 
maximum likelihood locations for the low- and high bands and 
the nominal position of r\ Car are <4igh-iow = 0.033°, dhigh-nominai 
= 0.038°, diow-nominai = 0.012°. The centroid of the low-energy 
band is well inside the 95.4% error region of the high-energy 
band. Although the high band's centroid lies slightly outside the 
low band's 99.7% error circle, the error regions substantially 
overlap. The difference to the indicated catalogue position of 
2FGL J1045.0-5941 is understood by acknowledging the dif- 
ference in exposure and the decomposition into low- and high- 
energy bands. We stress that all aforementioned positions suffi- 
ciently agree to allow interpretation in a single- source scenario. 

Considering the possibility of observing the y-ray emission 
superimposed by two sources in close spatial coincidence, we 
consider this to be extremely unlikely. Although the nominal 
distance between the centroid of the observed low- and the high- 
energy emission would not exclude source confusion between a 
low-energy source with exponential cutoff and a hard- spectrum 
source at high energies, the probability of finding two spatially 
coincident sources that exhibit a similar flux modulation pattern 
(see Section [33]), characteristic only for objects in the rare class 
of long-period y-ray binaries, appears to be the main argument 
for rejecting this hypothesis. 



3.2. Spectral analysis 

The energy spectrum of 2FGL J 1045.0- 5941 has been deter- 
mined by a maximum-likelihood fitting technique using the tool 
gtlike (dataset and source model as described above). The energy 
bins were chosen such that a minimum detection significance of 
10 cr per bin is required (logarithmic scale with three bins be- 
tween 0.2 GeV and 1 GeV and five bins between 1 GeV and 10 
GeV). Owing to statistics, the criterion was lowered to require at 
least 5 cr detections per energy bin in the range from 10 GeV to 
100 GeV. The highest energy bin (100 GeV to 300 GeV) corre- 
sponds to a 3 cr detection significance. For each bin a single PL 
shape was assumed. The normalization parameter and the index 
were allowed to vary during the fitting process. The effective en- 
ergy of each bin is defined as the midpoint of the integral photon 
flux in the bin; therefore each energy bin carries some informa- 
tion about the underlying photon spectrum. 

A spectral fit over the energy range (0.2 to 100 GeV) was ob- 
tained by fitting a source model with a CPL+PL function (both at 
the catalogue position of 2FGL J1045.0-5941) to the dataset up 
to 100 GeV using the maximum-likelihood tool gtlike. Events at 
energies E > 100 GeV have not been considered in fitting spec- 
tral models, because they do not significantly influence the re- 
sults. Figure [2] shows the energy spectra of 2FGL J1045. 0-5941 
along with the overlaid broad-band spectrum. The CPL+PL 
m ulti-compon ent functional form has been applied previously 
by lAbdo et al.l {2010b) and continues to yield better spectral rep- 
resentation than a single PL or CPL function. For comparison to 
a signal without a second emission component, the spectrum of 
the nearby pulsar PSR J1048-5832 (1.2° angular distance) is 
shown, too. The best- fit parameters are a spectral index T=1.97 
± 0.05 and energy cutoff E cuto ff=3.5 ±0.1 GeV with an inte- 
grated flux of F^ooMev =(1 20 + 0.04) x 10" 7 cmV 1 for the 
CPL, which dominates the low-energy band and T=1.94 ± 0.04 
and an integrated flux of F^MeV =(3 2 + q.3) x 10 " 8 cmV 1 
for the PL dominant in the high-energy band. The correspond- 



ing formal significances of the two spectral components are 30 cr 
and 16 cr, respectively. 

3.3. Temporal analysis 

We now investigate the flux time history of the observed y-ray 
emission. To search for coarse temporal signatures, two adja- 
cent energy intervals were analysed: 0.2 to 10 GeV and 10 to 
300 GeV, respectively. This choice corresponds to the apparent 
structure in the spectral energy distribution (SED) at 10 GeV 
as evident in Figure [2 which marks the transition between two 
spectral components. On the basis of 35 months of Fermi-LAI 
data it is not useful to search the power spectrum for the ex- 
istence of the characteristic ~5.5 yr orbital period in the r] Car 
system. Therefore we restrict the temporal analysis to flux vari- 
ability investigations, conducted over the whole energy interval 
and in the two energy regimes separately. 

3.3.1 . Flux studies in the lower-energy band (0.2 to 1 GeV) 

In Figure [3J left, we show the flux time history of 
2FGL J1045.0-5941 for the energy band 0.2 to 10 GeV as ob- 
tained by likelihood analysis. The dataset was divided into four- 
teen consecutive intervals, each representing 2.5 months. For 
each time bin a CPL function (all parameters free) was fitted to 
the data. The average photon flux is (1.29+0.01) x 10" 7 cm'V 1 . 
Taking this value as the hypothesis for a non- varying source, a 
;^ 2 -test gives a probability of 99.986% (corresponding to 3.8 cr) 
that this hypothesis is false. Conducting the same study on the 
two nearby pulsars PSR J1048-5832 and PSR J1044-5737 (at 
angular separation of 1.2° and 2.0°), which are expected to show 
a steady flux behaviour, we obtain merely a 13 cr and 0.8 cr 
significance of null-hypothesis violation. 

We conclude that in contrast to the nearby pulsars, the y- 
ray source 2FGL J 1045.0 -5 941 does show a significant degree 
of variability in the low-energy band. Detailed analysis shows 
that this is mainly because of the low flux values obtained for 
the most recent 12.5 months of our data sample. The mean flux 
for this period lies about a factor 1.7 below the mean flux of 
the remaining data sample. The nearby pulsars do not show this 
downward trend. 

The corresponding cumulative TS -value evaluation for 
2FGL J1045. 0-5941 is shown in Figure[3j right. A non-variable 
source is expected to s how a linear increase in the cumulative 
TS-value ( Matt ox et al.l [T 996) and deviations from that can be 
an indication of statistical fluctuations, stochastic flux variabil- 
ity, or periodic flux variation. In addition to minor deviations, the 
figure indicates an approximately linear rise in the cumulative 
TS-value for the first 22.5 months of 2FGL J1045.0-5941 data 
(blue dotted line). Then, however, the growth in TS decreases in 
correspondence to the lower flux values in Figure |3] left. Thus, 
the conclusions motivated by the light curve are reinforced by 
the cumulative TS evaluation. 



3.3.2. Flux studies in the higher-energy band (10 to 300 
GeV) 

Figure |4] is the equivalent of Figure [3] for the energy regime be- 
tween 10 and 300 GeV. Due to the lower statistics at higher ener- 
gies, the time binning was increased from 2.5 to 5 months. The 
first two bins - representing the flux of the first ten months of 
Fermi-LAI observation - are about a factor of 2.5 higher than 
the following five bins. For the most recent five months only 
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an upper limit could be obtained. A ^ 2 -test shows that a null- 
hypothesis of observing a constant flux can be rejected with 
7.4 <t significance in favour of flux variability. 

This conclusion is supported by the evaluation of the cor- 
responding cumulative TS -values. Due to the cut-off spectra of 
the nearby pulsars there is no other source in the vicinity of 
2FGL J1045. 0-5941 that can serve for comparison in this en- 
ergy interval. Yet, the behaviour of the TS -value is indicative 
of y-ray variability: whereas the cumulative TS- value constantly 
rises seemingly linearly and steadily over the first ten months 
(blue dotted line), it flattens and turns over afterwards. A line 
connecting the origin of the graph and the last cumulative TS- 
value (green dashed line) would indicate steadily increasing de- 
tection significance over time and thus absence of flux variabil- 
ity. Because this is not the case, we have to conclude that the 
cumulative TS -curve is indicative of flux variability. However, 
this indication builds up only after the first ten months of data; 
after an initial period of constant flux, the y-ray emission at en- 
ergies 10 to 300 GeV begins to fade out. 



3.3.3. Spectral evolution 

Figure [5] shows two SEDs (determined as in Section IX2t . The 
SED on the left represents the first ten months of our dataset 
for which an increased flux was observed at energies 10 to 300 
GeV. The SED on the right covers the remaining 25 months. As 
suggested by Figureffl the data points of the hard component ap- 
pear to be systematically lower in the latter part of the data set. 
Similarly as in Section IX2l a CPL+PL was fitted to the data for 
both time intervals using energies below 100 GeV. The most sub- 
stantial change is to be seen in the high-energy PL component 
in which the index drops from T=1.41 ± 0.08 for the first ten 
months to T= 2.06 ± 0.06 for the latter 25 months. Using the ob- 
tained spectral parameters to derive the integrated flux from 10 
to 100 GeV (the contribution of both spectral functions are taken 
into account), we obtain F^jggS^v = (1.2+g-|) x 10" 9 cmV 1 and 

F <ioogIv = (°- 49 ± °- 05 ) x 10 ' &9 cm " 2s-1 > which implies a mod- 
ulation by a factor of about 2.4. The difference in flux therefore 
agrees well with the findings of Section IT. 3. 21 



3.3.4. Connection to the orbit of the rj Car system 

As indicated by the vertical line in Figure [4] the period of in- 
creased flux in the high-energy component (10 to 300 GeV) at 
the start of Fermi-LAI observations corresponds to the perias- 
tron passage of the n Car binary system. To further illustrate this 
relation, Figure^ left, shows the flux of 2FGL J1045.0-5941 at 
energies from 10 to 300 GeV over the orbital phase of n Car. The 
periods of increased flux clearly correspond to orbital phases 
close to periastron. Figure right, shows flux in relation to the 
physical distance of the two stellar components (ranging from 
1.5 AU at periastron to -30 AU at apastron, which was reached 
by October 2011). The intervals of increased flux correspond to 
a time in which the separation of the stars is less than about 12 
AU, with a flux maximum (bin 1) just prior to the periastron 
passage. Otherwise the data indicate a clear decrease in flux as 
the binary components approach apastron. The same is shown in 
Figure |7] for the low-energy component (0.2 to 10 GeV). In this 
case the flux decrease seems to commence later, with periods of 
lower flux corresponding to stellar separations of more than 25 
AU. 



4. Discussion and summary 

As the analysis above has shown, we see a clear flux decrease in 
relation to the orbital state of n Car. The decrease is significant 
at low- (0.2 to 10 GeV) as well as high- (10 to 300 GeV) en- 
ergies, albeit showing a different time dependence (see above). 
In general, the observed flux decrease can be qualitatively un- 
derstood with predictions in the framework of colliding-wind bi- 
nary models for high-energy y-ray emission (Reimer et al. 2006; 
iBednarek & Pabichll201 ll) . As the stars move away from each 
other, the matter density - as well as the radiation density - de- 
creases in the wind collision region where particle acceleration 
and subsequent y-ray emission are thought to occur. Although 
the observed flux variation does not yet allow an unambiguous 
association, we consider in the following that the detected emis- 
sion indeed arises from n Car for reasons explained above. 

Note that source modulations of the two spectral components 
at the same amplitude can be fully recovered for the two energy 
intervals only if the model used for the diffuse gamma-ray emis- 
sion and for nearby point sources is perfect (e.g., without statisti- 
cal and systematic uncertainties). Any realistic assessment of the 
backgrounds (e.g., for statistical fluctuation or, even more, sys- 
tematic deviations between model and reality) will have an im- 
pact on the principal ability to recover source modulation. This 
will matter comparably less at higher energies due to the better 
signal-to-noise ratio. Hence, our inability to produce a perfect 
diffuse emission model principally limits our capability to re- 
cover variability in the low-energy component of n Car. 

The presence of phase-lo cked flux changes in any emission 
component refutes the idea of Ohm et al. (2010), who suggested 
that a y-ray signal from r\ Car can be related to particle accel- 
eration in the expanding blast wave that originated in the great 
eruption of 1843 when the primary star (the LBV) ejected about 
10% of its total mass in a massive outburst. Through shock inter- 
action with the surrounding interstellar medium the blast wave 
that was created at that event might produce charged particles 
and subsequent y-ray emission. In such an emission scenario one 
would not expect to see y-ray flux variability in connection with 
the orbital configuration of the n Car binary system. 

Based on the observations with the LAT of two appar- 
ently distinct emissi on components, a model has been proposed 
(Farni eret al.l 1201 lb to identify the low-energy component as 
due to inverse Compton scattering of n Car's photospheric pho- 
ton field off a relativistic electron population extending up to 
~10 4 MeV, and the high-energy component above ~10 GeV as 
due to hadronic interactions of a relativistic proton component 
extending up to ~10 4 GeV with a high-density material in the 
shock region (enhanced by a factor ~10 with respect to the un- 
perturbed wind density). This requires electron and proton ac- 
celeration at maximum rate (Bohm diffusion regime) for values 
of the shock magnetic field in the sub-Gauss range. The target 
material and photospheric photon density at the shock location 
change with orbital phase as the inverse square of the distance 
between the shock and the LBV. Therefore, phase-locked orbital 
variations of both inverse Compton and hadronic components 
are expected unless the intensities of the putative accelerated par- 
ticle populations change along the orbit such that it counteracts 
this trend. We stress here that both components are expected to 
vary significantly with the same amplitude for an isotropic target 
radiation field and no significant orbital variations in the accel- 
erated electron-to-proton number densities. If anisotropy effects 
were considered, even larger variability amplitudes are expected 
in the inverse Compton component than in the hadronic com- 
ponent. Applied to the energy spectra determined for the first 
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10 months and the subsequent 25 months of LAT observations 
(see Figure O for which the mean stellar separation were ~8 AU 
and -21.8 AU, respectively, the flux amplitude is expected to 
decrease from the first to the second time period by more than a 
factor 7 in both spectral components, possibly by an even greater 
factor for the leptonic component. The observed spectra, how- 
ever, reveal flux variations in the high-energy component by a 
factor 2.5 while the low-energy component shows a rapid spec- 
tral decline at a few GeV with no indication for a phase depen- 
dence of the spectral shape, spectral index, or cutoff energy. 

We note that the observed low variability of the low- 
energy y-ray component with phase-locked variations of the 
high-energy component can be plausibly interpreted as inverse 
Compton scattering in the photospheric radiation field of n Car 
and hadronic interactions with the wind material, respectively, if 
the accelerated electron number density is higher around apas- 
tron than periastron while at the same time the accelerated proton 
number density does not exhibit significant orbital variations. 

Motivated by the peculiar shape of the measured spectra and 
the complex strong radiation fields in the vicinity of n Car, we 
here consider yet another scenario. Photon-photon absorption 
has long been suspected to play a significant role in the dense 
radiation environments of massive stars. For instance, first in- 
dications of y-ray absorption modulating a Te V spectrum have 
been reported from the y-ray binary LS 5039 (lAharonian et al.l 
l2005h . 

For n Car, the observed y-ray spectra can be well represented 
by a primary differential photon spectrum oc E~ 2 exp(-E/E cut ) 
with E cut ~ 250 - 500 Ge V that suffers from y-ray a bsorption in 
a black-body photon field (iGould & Schrederll 19671) of effective 
temperature ~ 5 • 10 5 K (corresponding to target photon ener- 
gies of ~ 116 eV at the maximum of the black-body radiation 
field) during the first ten months of observations, and decreases 
to ~ 2 • 10 5 K (~ 46 eV) during the subsequent 25 months (see 
Figure [8]). The required maximum optical depths are r yy ^ max « 



1 . 1 and r y 



2, corresponding to fluxes (at the distance of 



n Car) of ~ 6 • 10" 8 erg cm~ 2 s _1 and ~ 9 • 10" 11 erg cm" 2 s _1 , 
respectively. Considering the total spectrum (35 months), the 
required black-body absorber has an effective temperature of 
~ 3.1 • 10 5 K (~ 72 eV) and r yy , max « 1.3 corresponding to a 
flux of ~ 2 • 10" 9 erg cm~ 2 s _1 . In this scenario, a possible ori- 
gin of the absorber could be hot X-ray gas, as observed recentl 
by XMM-Newton an d Chandra (e.g., iHamaguchi 



/ea re cently 
et al.l I2007J: 



iTownsley et al.ll20lll) . Both experiments revealed spatially ex- 
tended, structured X-ray emission components that surround the 
binary system. Figure (left) illustrates the geometry of such 
an external absorber scenario. The change of the LAT spectrum 
would indicate a slightly hotter rarefied plasma towards perias- 
tron than towards apastron. 

In a second scenario the absorber could be linked to n Car's 
wind collision (e.g., hot shocked gas). In this case (see Figure[9] 
(right)) the angular dependence of the absorption coefficient and 
threshold energy on the photons' collision angle has to be taken 
into account. The resulting optical depth therefore depends on 
the orientation of the wind collision region, the line of sight, 
and the stars of the n Car binary system. For example, for a 
mean angular phase for the first 10 and subsequent 25 months 
of 0.973 and 0.422, respectively, and an inclination angle of 45° 
and to = 285°, the propagation angle (between the collision re- 
gion and line of sight) for target photons impinging from the 
direction of n Car is ~ 73° and ~ 119°, respectively. In this case 
we find satisfactory spectral representations for an absorber of 



temperature T ~ 10 6 K and r r 



1 for the first 10 (see 



Figure [U left), and T ~ 2.5 • 10 5 K and r yy , max « 2 for the latter 
25 months (see Figure [8] right). 

Orbitally modulated y-ray emission up to TeV energies has 
been observe d from a numbe r of high-mass X-ray binary sys- 
tems (see e.g. jHill et al. 2010) with orbital periods ranging from 
3.4 years (PSR B 1259-63) down to a few weeks or days (e.g., 
LSI +61 303, LS 5039). For the latter two systems, peak fluxes 
at MeV/GeV energies occur close to superior conjunction where 
the companion is behind the star (which also turns out to be close 
to periastron passage), while the flux decreases towards inferior 
conjunction/apastron. This behaviour is similar to the observed 
flux modulation in the spectrum of n Car despite the lack of de- 
tection at VHE y rays. 

Ultimately, if one concludes that the observed decrease in 
y-ray flux indeed relates to orbital conditions in the n Carinae 
system, a prediction can be made for how the y-ray flux will 
evolve until periastron is reached again. It is expected to remain 
at a low level until the end of 2013, when the flux should increase 
again and, ultimately, reach maximum shortly before the next 
periastron passage in mid-2014. The observation of this regular 
variability pattern will be a decisive test for current models of 
y-ray emission in colliding-wind binary systems. Until then, we 
anticipate that numerical simulations will consider the multitude 
of details not yet taken into account in present y-ray emission 
models (e.g., complexity in stellar radiation fields specific for 
the 7] Car system and the Homunculus nebula). This will allow 
a quantitative comparison of model predictions to the observed 
y-ray data over a full orbit. 
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Figure 1. Test statistic (TS) maps of the celestial region around 2FGL J1045. 0-5941 as obtained by likelihood analysis. The TS- 
value corresponds to the likelihood of the source being located in a specific grid point. Left plot: 0.2 to 10 GeV, right plot: 10 to 
300 GeV. The confidence contours (green for the high band, white for the low band) mark the 68.3%, 95.4% and 99.7% uncertainty 
regions for the location of the source of interest. The dashed yellow circle marks the catalogued 95% source location region of 
2FGL J1045. 0-5941. The nominal position of r] Car is indicated by the yellow cross. 
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Figure 2. Spectral energy distributions of 2FGL J1045.0-5941 (red dashed) and PSR J1048-5832 (green solid) obtained by 
likelihood analysis. The spectral model is obtained by fitting a source model with CPL and PL components at the location of 
2FGL J1045. 0-5941 to the dataset (35 months, 0.2 < E < 100 GeV). The error bars are oiler type. The upper limits (represented 
by arrows) were determined such that the difference of the logarithmic likelihood values (with and without an additional trial point 
source at the indicated flux value) corresponds to 1 <x. 
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Figure 3. Left: Flux time history of 2FGL J1045. 0-5941 in the 0.2 to 10 GeV energy band obtained by likelihood analysis. Each 
bin represents 2.5 months of data. The fl ux error bars are o f 1 cr type. The dashed vertical line indicates the time of the periastron 
passage of the // Car binary system (e.g., Parkin et al. 20091). 

Right: Cumulative TS-value evaluation of 2FGL J1045. 0-5941 as obtained by likelihood analysis in the 0.2 to 10 GeV energy band. 
Each data point represents the cumulative TS-value of the source for the time interval from the start of the Fermi-LKY data taking. 
The dotted blue line is a linear fit to the first nine data points. The dashed green line connects the origin and the cumulative TS-value 
obtained last. 
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Figure 4. Left: Flux time history of 2FGL J1045. 0-5941 obtained by likelihood analysis. Each bin represents five months of data 
in the 10 to 300 GeV energy band. The flux error bars are of 1 cr type. The upper limit was determined as described in Figure The 
dashed vertical line indicates the time of the periastron passage of the r] Car binary system. 

Right: Cumulative TS-value evaluation of 2FGL J1045. 0-5941 as obtained by likelihood analysis in the 10 to 300 GeV energy 
band. Each data point represents the cumulative TS-value of the source for the time interval from the start of the Fermi-LKT data 
taking. The dotted blue line is a linear fit to the first four data points. The dashed green line connects the origin and the cumulative 
TS-value obtained last. 
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Figure 5. Spectral energy distributions of 2FGL J1045. 0-5941 for the first 10 months (left) and latter 25 months (right) of the 
dataset obtained by performing likelihood analysis. The spectral model is the sum of CPL and PL components at the position of 
2FGL J1045.0-5941 to the 0.2-100 GeV data. 
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Figure 6. Flux time history of 2FGL J1045. 0-5941 as determined by likelihood analysis for ranges of the orbital phase of r/ Car 
(left) and the physical separation of the stellar components (right) for energies 10 to 300 GeV. Each data point represents five months 
of data. Because periastron is traversed very quickly, the data points around orbital phase are very broad. The region between the 
two dashes at the left of the horizontal error bar of the second data point on the right graph is traversed twice. 
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Figure 7. Flux time history of 2FGL J1045. 0-5941 as determined by likelihood analysis for ranges of the orbital phase of ?] Car 
(left) and the physical separation of the stellar components (right) for energies 0.2 to 10 GeV. Each data point represents 2.5 months 
of data. Because periastron is traversed very quickly, the data points around orbital phase are very broad. The region between the 
two dashes at the left of the horizontal error bar of the leftmost data point on the right graph is traversed twice. 
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Figure 8. Representation of SEDs (as in Figure overlaid with the unabsorbed photon spectrum oc E~ 2 exp(-E/E cut ) with E cut ~ 
250 - 500 GeV (thick line) and one that suffers from y-ray absorption (dashed line) shown for the first 10 months (left) and the 
subsequent 25 months (right). 
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Figure 9. Illustration of the projected orbit of the binary system considering y-ray emission towards the line of sight and its possible 
attenuation by y-ray absorption as described in the text. Scenario 1 (left): external black-body absorber (e.g., hot gas surrounding 
the binary system). Scenario 2 (right): internal absorber (e.g., hot shocked gas in the wind collision region.) 
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